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CHAPTER I 



INTRODUCTION 

Business men the world over, have long recognized 
that records and reports of the physical quantities and 
costs of inventories, including their balances, changes and 
relationships, are essential to managerial appraisal of past 
business performance and to planning and control of future 
operations.^ The objective of inventory management is to 
maintain the investment in inventories at the lowest amount 
which is sufficient to meet production, sales and financial 
requirements of the enterprise. The inventories must be 
adequate to maintain an efficient level of operations and to 
meet, within reason, the needs of customers. ' However, it 
must not be greater than is necessary to meet these require- 
ments because of interest costs, the cost of handling and 
storing excessive quantities of inventory, the dangers of 

adverse price changes and obsolescense, and the increased 

2 

exposure to physical deterioration. Thus the requirement 
for an accurate inventory determination is a prerequisite to 
effective inventory management. 

For most business enterprises, the determination of 
inventory assets presents no unusual problems since most 
inventory assets can be physically ascertained. However, 
in certain businesses such as the manufacture of chemicals, 
major portions of their inventory assets are in bulk raw and 
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finished materials. The problem of accurately determining 
inventory balances of bulk materials has continually plagued 
business managers because such materials are not susceptable 
to recognized mensuration techniques. The quantifying of 
bulk inventories is especially difficult when the economics 
of handling and storing the bulk materials dictate that they 
be stored in stockpiles which have a great degree of surface 

3 

non-conformity and no constraining boundaries. 

Purpose and Scope 

The purpose of this thesis is to determine, by experi- 
ment, the minimum measurements necessary to estimate the con- 
tent of irregular solids within a prescribed error using nu- 
merical analysis methods. The work encompasses a laboratory 
study utilizing granular ammonium nitrate for construction of 
stockpiles having varying configurations. The study and the 
results of the laboratory tests are to be applied to determine 
the reserves of ammonium nitrate stored in the bulk storage 
warehouses of the Cooperative Farm Chemical Association at 
Lawrence, Kansas. 

In the laboratory experimental project, a digital 
computer program and a manual computation method are de- 
veloped to estimate the content of stockpiles of granular 
materials. Additionally, the angle of repose of coated 
ammonium nitrate granules is ascertained and compared with 
the accepted value for this characteristic. 
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Historical Perspective of Volume Calculations 

For centuries, mathematicians have known the formulas 

4 

for calculating the volume of linear and non-linear solids. 
These formulas require explicit measurements of the length 
of a line or the size of an angle. Examples of such volumes 
are spheres, cones, pyramids and solids of revolution such 
as ellipsoids and cylinders. Volume determinations may also 
be based on implicit measurements. This field encompasses 
the science of integral calculus wherein volumes are computed 
by integration of either a known analytic function or a fin- 
ite numerical difference expression of the area or one linear 
dimension as a function of the other two dimensions. However, 
very little application has been made using the numerical 
difference approach since the amount of. computational work is 
large by this method. Thus, until the advent of the digital 
computer, the best course of action was to estimate the volume 
of irregular solids using the classical explicit volume form- 
ulas and to accept the error inherent with this method. The 
present study is an investigation of the benefits to be 
gained through the use of numerical integration of various 
numerical difference expressions of the area or height of a 
stockpile. 

Problems .Associated with Stockpile Volume Estimation 

Undoubtably the biggest problem in measuring stock- 
piles is their size. This characteristic coupled with the 
irregular shape of the stockpile surface precludes normal 
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measurement methods based on geometrical formulas- Addition- 
ally, if the stockpile is composed of granular material it is 
impractical, if not impossible, to take depth measurements 
over the surface because of the fluidity of the granules. 
Consequently, the practice of surveying the surface with a 
transit came into vogue in an attempt to estimate the content 
of the stockpile. This procedure is quite laborious and the 
accuracy of the volume calculations leaves something to be 
desired . 

Aerial photocontouring is feasible for outdoor stock- 
piles. Results by this procedure are fairly accurate but 
this method is expensive and takes considerable time due to 
the limited number of photocontouring firms. 

Photocontouring on The other hand is impractical for 
many small outdoor industrial stockpiles. Needless to say it 
is nearly impossible to photocontour indoor stockpiles due to 
factors such as size of the warehouse and environmental con- 
ditions such as dust. 

Another problem associated with determining the con- 
tent of stockpiles is the amount of time and cost involved 
in taking the measurements and performing the necessarily cal- 
culations to fix the inventory quantity. In exercising sound 
inventory management control, engineers and businessmen 
should be knowledgeable of the alternatives available in terms 
of inventory accuracy versus the time and costs required to 
attain that accuracy. 



CHAPTER II 



NUMERICAL ANALYSIS METHODS 



Numerical Integration 

Integration may be defined as a process of summation. 
Numerical integration is defined by J. B. Scarborough "as 
the process of computing the value of a definite integral 
from a set of numerical values of the integrand".^ This 
process is sometimes called mechanical quadrature when 
applied to the integration of a function of a single vari- 
able; when applied to the calculation of a double integral 
of a function of two independent variables it is called 
mechanical cubature . 

For numerical integration, the problem is to estimate 
the numerical value of the integral: 



I(x) 




f(x) 



dx 



( 2 . 1 ) 



when f(x) either is too complicated an analytical function 
to permit the integration in an analytic manner or when the 
value of f(x) is given only in tabulated form. One method 
of solving the problem is to expand f(x) in an infinite 
series, the individual terms of which can be integrated. 
Thus if 

OO 

f(x) 0 ^ (X) 

n=0 



( 2 . 2 ) 
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and one can determine the constants 



n 




then 



b 

I(x) = J f(x) dx = ^ (2.3) 

a A 

n=0 

Such a method is practical, provided this series converges 
rapidly, i.e., only a few terms need to be calculated to 
determine I(x) to any desired accuracy. 

The problem of numerical integration is solved for 
the general case by representing the integrand by a suitable 
polynomial, for a given interval and subsequently integrat- 
ing the polynomial between the desired limits. The poly- 
nomials used for this purpose are called interpolating poly- 
nomials. They are derived using the methods of finite dif- 
ferences. Finite difference methods are covered very exten- 
sively in the literature of numerical analysis. An excellent 
presentation is that by Scarborough. An important property 

of the usual power series interpolating polynomial is that 

"t h "t h 

the n^ differences of a polynomial of the n^ order are con- 
stant when the values of the independent variables are taken 
in arithmetic progression, that is, at equal intervals. This 
property makes it possible to represent any function, includ- 
ing one in tabular form, by a polynomial if its differences 
at some degree become constant or nearly constant. 



When the tabulation is made at equal intervals, the 
accuracy of the quadrature formula depends upon the size of 
the interval, the limits of integration, and the degree of 
the interpolation formula. 



General Quadrature Formula 

Scarborough and others present a derivation of Newton's 
formula for forward interpolation which is the basis for de- 
veloping a general quadrature formula for equidistant ordi- 

. 7,8,9 

nates. ’ ’ 

This general quadrature formula is: 

(2.4) 



x^+nh 



3 2 



/"O ““ r 2 3 2 

J ydx = h ny^ ^ Ay^ + 

X L. 
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+ (-7 n +n ) — + (.-= 



3._2^ '^•'^0 . 3n' 



3i 



5 2 



lln^ ^2. ^ ^0 
-3 3n ) 



. ^n- ^.5 . 35n"^ 50n^ . 2^ ^ 

+ (-^ - 2 n + j- + 12 n ) — jj- 



, .rJ 15n^ 5 225n"^ , 2740^ 2. ^^^0 






where n is the number of equidistant intervals of width h, 
and A is the difference operator on the value of the vari- 
able y. 

The notation used for the difference operator in the 
general quadrature formula can be explained as follows. If 
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Yq, y^, 72 ••• Yj^ denote a set of values of any function 

y = f(x), then y^ - y^, - y^, 73 - 72 ’ •** Vn “ '^nl’ 

called the first differences of the variable y. Denoting 
these differences by ^7 q> have ^7 q ~ 

" ^0’ ^^1 = ^2 " >^1’ ^^n-1 = ■ ^n-l’ = ^n^-l ” "^n* 

The differences of the first differences are called the 

2 2 

second differences and are denoted by A Yq, A y^, etc. 

A^Yq = Ay^ - Ayq = 72 - 2y^ - y^ 

A^Yi = AY 2 - Ay^ = 73 " - y^ 

In the like manner, the third differences are; 

A^Yq = A^Yi - = Y 3 - 372 " ^^1 " ^0 

A^y^ = A^Y2 - A^y^ = Y4 ~ 373 - 872 - y^ 

Numerical analysis literature is replete with diagonal 
and horizontal difference tables. The use of these tables 
depends on which method of interpolation is to be used. In 
this study, forward interpolation will be used since the 
tabulated data values begin and end with a zero value over 
the range of integration. 

From the general quadrature formula, specific quad- 
rature formulas may be obtained by varying the value of n. 

Three quadrature formulas have found widespread use 
due to their simplicity and accuracy. These formulas are; 
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Trapezoidal Rule, Simpson's Rule and Weddle’s Rule. They 
are obtained from equation (2.4) by letting n = 1, 2, and 6 
respective ly . 



Trapezoidal Rule 

Substituting n = 1 in equation (2.4) and neglecting 
all differences above the first, one obtains 



The geometric significance of this formula is that a given 
function y(x) is replaced by n/l arcs (straight-line segments), 
i.e., first degree polynomials. It is exact if y = f(x) is a 
polynomial of degree 1. 

Simpson’s Rule 

Substituting n = 2 in equation (2.4) and neglecting 
all differences above the second, one obtains 



Simpson’s Rule corresponds to the representation of the given 
function y(x) by n/2 arcs of a second degree polynomial. It 
is exact if y = f(x) is a polynomial of degree 2 or less. 
While it is not readily apparent, Simpson’s Rule is exact if 
y = f(x) is a polynomial of degree 3 because the third dif - 
ference ( A y^) is zero and the first term dropped involves 




0 



(2.5) 







( 2 . 6 ) 
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Weddle's Rule 

Substituting n = 6 in equation (2.4) and neglecting 
all differences above the sixth, one obtains 
»x +nh ^ 

^0 



(2.7) 

The geometric interpretation of this rule is that a given 
function y(x) is replaced by n/6 arcs of a fifth degree poly- 
nomial. This rule is exact for fifth-degree polynomial or 
lower. This formula requires at least seven consecutive 
values of the function, or, stated alternatively, the inte- 
gration interval must consist of six subintervals of equal 
width h. 



Volume Computations 

To calculate the volume of any solid using quadrature 
formulas, it is only necessary to establish an x-y grid sys- 
tem and to determine the values of the Z coordinate corre- 
sponding to each point of that system. Figure (2.1) shows 
the graphical representation for this method of mechanical 
cubature. Thus, if a solid is divided into cross sections, 
each ordinate of that cross section may be measured and tabu- 
lated. These values may then be used in any acceptable quad- 
rature formula such as Simpson’s or Weddle’s formulas to 

determine the area, A 

’ cs 



of each cross section. 
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Figure 2.1 

Graphical Representation of 
Mechanical Cubature on a Irregular Solid 




Volume = 



Acs dx 



I. 
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Ac3 = J Z dy 



( 2 . 8 ) 



Subsequently, the cross sectional areas, may be used as 
the ordinates of the length function and themselves be inte- 
grated using any acceptable quadrature formula. 



In this manner it is possible to compute the volume of a solid 
nine ways if one is limited to the three most frequently used 
quadrature formulas based on equal spacing of the independent 
variable. Generally, the accuracy of the double quadrature 
method depends on the number of data points and on the shape 
of the surface of the solid. An important characteristic of 
all the quadrature formulas based upon an interpolating poly- 
nomial is that the formulas are numerically exact for those 
instances in which an n^ order polynomial y(x) is used and 
for which instance the actual function being approximated is 
of order n or less. 



Volume 




( 2 . 9 ) 



c 
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CHAPTER III 



EXPERIMENTAL PROCEDURE 



Apparatus 

The experimental apparatus shown in Figure (3.1) was 
utilized to construct and measure the different configura- 
tions of stockpiles. 

The test site consisted of a 4' x 8’ sheet of 3/4” 
plywood mounted on a large stable table. One inch ruled co- 
ordinate paper was then attached securely to the surface of 
the test bed. A rectangular area 24” x 60” was ruled on the 
coordinated paper to establish the boundary of the experi- 
mental stockpiles. These boundaries were used as the ref- 
erence lines in measuring the size of the stockpile at its 
\ 

base . 

Conveyor System 

An overhead track system was devised to simulate an 
overhead conveyor system. A small trolley, carrying a steel 
funnel, was allowed controlled movement along the overhead 
track. The distance from the funnel spout- to the test bed 
was designed to be eighteen (18) inches. The model conveyor 
duplicated the formation of stockpiles as they appear in the 
plant warehouse of the Cooperative Farm Chemical Association. 

Measuring Device — Depth Gauge 

Two laboratory point-depth gauges. Lory Type - A, manu- 
factured by Leupolo and Stevens Instruments, Portland, Oregon 
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Figure 3.1 

Experimental Apparatus for 
Constructing Laboratory Scale Stockpiles. 







1 
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and capable of measuring distances to an accuracy of .001 
foot, were mounted in two ”T" frames which were constructed 
from two inch angle iron. This mechanism permitted the 
measurement of the stockpile surface at one inch intervals 
both across and along the length of the stockpile. The 
gauges were leveled and calibrated for each series of cross 
sectional measurements. 

Materials 

The ammonium nitrate used was taken from the standard 
80 lb. bags supplied by the Cooperative Farm Chemical Associ- 
ation. A standard 80 lb. bag of ammonium nitrate was de- 
termined by repeated measurements to contain an average of 
2929.0 cubic inches. The initial density was determined to 
be 47.19 Ibs./cu.ft. Frequent checks were made to determine 
if the stockpile compacted during the four weeks duration of 
the experiment. These measurements indicated that possibly 
a very slight compaction occurred mostly near the center re- 
gions of the stockpile. However the change in height of the 
pile at all positions was within the experimental error in- 
herent in the height measurements. Therefore, it was con- 
cluded that the density change was not sufficient to warrant 
its inclusion in the volume calculations. The original volume 
was therefore used as the basis for all calculations during 
the experiment. 

Miscellaneous Equipment 

A stone jar was used to measure the volume of the 
ammonium nitrate used in the investigation. The volume of 
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the jar was determined by measuring the volume of the water 
necessary to fill it. A calibrated 1000 cc cylinder was used 
to fill the jar. The capacity of the stone jar was deter- 
mined to be 842.0 cu.in. 

Measurement Procedure 

A preliminary series of measurements was conducted 
using a single small conical pile of ammonium nitrate granules, 
to test the proposed measuring procedure. Initially, it was 
proposed that the tip of the depth gauge would contact a 
granule on the pile surface and thereby determine the height 
of the pile at that point. Measurements on the small pile 
using this procedure proved difficult and gave inaccurate re- 
sults. This inaccuracy was caused by the irregular arrange- 
ment of the ammonium nitrate pellets or granules. The tip 
would displace a granule and leave the tip in the air above 
the pile. After repeated measuring attempts, a technique was 
developed whereby the sharp needle-like point of the measure- 
ment probe penetrated the stockpile surface 1/16 of an inch 
(approximately one pellet diameter). This method resulted 
in consistent height measurements. 

The experiment began with the filling of the funnel 
with the ammonium nitrate granules. The nitrate was then 
allowed to drop continuously from the funnel conveyor in the 
overhead track system, falling the eighteen inches to the 
test bed. When sufficient granules had been deposited, such 
that the base of each sub-stockpile reached the ruled 
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boundaries, the funnel conveyor was then moved twelve inches 
to a new position. This procedure was repeated until the 
80 lbs. (2929 cubic inches) of ammonium nitrate was deposited 
on the test site. 

Figure (3.2) shows the configuration of the "undis- 
turbed” stockpile. The four concial shaped sub-piles were 
formed to approximate the typical storage situation in the 
warehouse. The conveyor system was then removed from the 
test site and the height measuring mechanism installed at one 
end. After calibrating the height gauge, measurements were 
taken over the surface of the stockpile at one inch intervals 
in compliance with the requirements of the quadrature formulas. 
In addition to the 1525 height measurements, 61 width measure- 
ments were taken of the outline which the base of the stock- 
pile made on the ruled paper. The angle of repose was also 
measured at various cross sections around the stockpile us- 
ing a jointed protractor. Table (3.1) tabulates these measure- 
ments . 

Upon completion of the measurements on the "undisturbed” 
stockpile, a second stockpile was created by selectively re- 
moving some material from the above mentioned "undisturbed" 
pile. It was felt that this would represent the configura- 
tion of a "working" or "disturbed" stockpile. A child’s toy 
shovel was used to remove small amounts of nitrate from several 
positions at the base of the undisturbed pile. This process 
simulated quite well the random removal of the ammonium ni- 
trate by a tractor loader. The amount of 842 cubic inches of 
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Figure 3.2 

Undisturbed Laboratory Scale Stockpile. 
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TABLE 3.1 

ANGLE OF REPOSE AT PEAK 
CROSS SECTIONS AS MEASURED BY PROTRACTOR 



Cross Section 


Undisturbed 

Stockpile 


Disturbed 

Stockpile 


Left 

Side 


Right 

Side 


Left 

Side 


Right 

Side 


13 


31.0 


31.0 


27.0 


28.5 


25 


31.0 


31.5 


28.5 


29.0 


37 


31.5 


30.5 


26.5 


30.0 


49 


30.5 


31.0 


28.5 


29.0 


END 


Left 


Right 


Left 


Right 




End 


End 


End 


End 




31.0 


30.5 


27.0 


27.5 


Average Angle 










of Repose 


30.9 


28.2 
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ammonium nitrate was removed from the base around the entire 
periphery of the original stockpile. Figure (3.3) shows the 
’’disturbed” stockpile. After the removal of the 842 cubic 
inches of material no attempt was made to alter the shape of 
the pile. As can be seen from Figure (3.4), the surface 
tends to be linear and the peaks have regressed into ridges. 
The ’’disturbed” stockpile was then measured in the same man- 
ner as was the ’’undisturbed” stockpile. 

Computational Models 

An explanation of the various computational models 
employed to ascertain the volume of the ’’undisturbed” and of 
the ’’disturbed” stockpiles will now be given. These models 
range in complexity from the siraplist in which the stockpile 
cross sections are assumed to be isosceles triangles to the 
most complex model wherein the depth measurements on a one 
inch square grid are numerically integrated using Weddle’s 
Rule. The computer programs for each model were devised to 
allow parameters such as angle of repose, and the number of 
measurements in the x and y directions of the x-y base grid 
to be considered as variables, to be read in as data. This 
flexibility permitted a study of the influence of these 
parameters on the volume calculation. 

The computer programs for the computational models 
used in this study may be found in Appendix A. 

Model # 1 

A computer program written in Fortran IV language 
was developed to compute the volume of stockpiles assuming 
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Figure 3.3 

Disturbed Laboratory Scale Stockpile. 
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Figure 3.4 

Disturbed Laboratory Scale Stockpile 



Depicting the Change in the Angle of Repose. 



% 
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isosceles triangular cross sections. The program calculates 
the cross sectional areas utilizing the angle of repose, a, 
as a parameter. Additionally, in computing the volume by 
numerically integrating the cross sectional areas over the 
length of the stockpile, three different quadrature formulas 
were used. The number of cross sections was also varied to 
observe its effect on the volume calculations. 

The equation for the area of the isosceles triangular 
cross section is shown below for the circumstance in which 
the angle of repose, a, and the base b are known 




tan a 

Cross Section Area 



h/| 

1/2 bh 

(l/2)(b)(^ tan a) 

„ 

4 



tan a 



( 3 . 1 ) 
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Model ft 2 

A computer program was developed to estimate the volume 
of a stockpile by numerical integration in two directions. 

The program calculates three different areas at each cross 
section utilizing the three quadrature rules: 

1) Trapezoidal Rule, 

2) Simpson’s Rule, 

3) Weddle's Rule. 

Also the number of depth measurements used in each cross sec- 
tional area can be designated externally and used as a model 
parameter. The program calculates the volume by numerical 
integration of the cross sectional areas previously obtained. 
For this phase of computation the number of cross sections 
and particular quadrature rule to be used can be selected 
externally and used as model parameters. 

Model # 3 

This computer program is a modification to Model #1 
and was written to specifically compute the volume of a stock- 
pile which is banked against a retaining wall. 

The cross sections are assumed to be isosceles triangles 
but with a portion of the triangle removed. The net area of 
each cross section is determined by computing the area of the 
isosceles triangle and then subtracting that part of the isos- 
celes triangle which projects beyond the retaining wall. 

Using the angle of repose as a parameter, the program 
proceeds to compute the volume of the stockpile by numerical 
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integration of the previously obtained cross sectional areas 
using the three quadrature rules. 

The equation for the area of a typical cross section 
is derived below for a banked stockpile. 



I 

I 

I 




Area(^^= |-(XLP)(Ph) 

= -|-(XLP)^ tan a 

Area(^^^^^= tan a 

Net Area = 2 x Area 1 - Area 2 

= (XLP)^ tan a - i(WF)^ tan a 



(3.2) 




CHAPTER IV 



ANALYSIS OF EXPERIMENTAL RESULTS 

General 

The results of the experimental study outlined in 
Chapter III will be presented in three parts. Part A is 
concerned with the computation of the volume and the deter- 
mination of the apparent angle of repose for the ’’undisturbed” 
stockpile. Part B describes the same results but for the 
’’disturbed” stockpile. Part C describes the results from the 
determination of the true angle of repose for coated ammonium 
nitrate granules. 

A. Undisturbed Stockpile Results 

Computational Models #1 and #2 were used in obtaining 
the volume and apparent angle of repose approximations. 

Model #1 

Table (4.1) tabulates the per cent error in the cal- 
culated volume obtained from computational Model #1 using 
as parameters; 

1. number of isosceles triangle cross sections 

2. quadrature rule 

3. apparent angle of repose. 

Because of the mathematical form of Model #1, the apparent 
angle of repose, a, is an independent parameter. The volume 
of the stockpile computed by the model is dependent upon the - 
value selected for a. If the volume of the ammonium nitrate 
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TABLE 4.1 



PER CENT ERROR IN CALCULATED VOLUME FOR 
UNDISTURBED STOCKPILE USING COMPUTATION MODEL ffl 



Cross 

Sections 


Per Cent Error in Calculated Volume 


Apparent 
angle of 
repose 


Trapezoidal 

Rule 


Simpson ' s 
Rule 


Weddle’ s 
Rule 




31.0 


5.49 


5.47 


5.45 




30.0 


1.36 


1.35 


1.32 




29.68. 


— 


— 


0.00 




29.67 


— 


0.00 


— 


61 


29.66 


0.00 


— 


— 




29.00 


-2.68 


-2.70 


-2.72 




28.0 


-6.65 


-6.67 


-6.69 




27.0 


-10.55 


-10.56 


-10.58 




31.0 


5.54 


5.20 


4.72 




30.0 


1.41 


1.08 


0.62 




29.85 


— 


— 


0.00 


31 


29.71 


— 


0.00 


— 




29.65 


0.00 


— 


— 




29.0 


-2.63 


-2.95 


-3.39 




28.0 


-6.60 


-6.91 


-7.33 




27.0 


-10.50 


-10.79 


-11.20 




31.0 


5.56 


7.51 


7.41 




30.0 


1.43 


3.30 


3.21 




29.65 


0.00 


— 


— 


13 


29.22 


— 


— 


0.00 




29.20 


— 


0.00 


— — 




29.0 


-2 . 62 


-0.82 


-0.91 




28.0 


-6.59 


-4.87 


-4.95 




27.0 


-10.49 


-8.83 


-8.92 




31.05 


0.00 








31.0 


-0.28 


6.28 


6 . 60 




30.0 


-4.18 


2.12 


2.43 


7 


29.48 


— 


0.00 


— 




29.41 


— 


— 


0.00 




29.0 


-8.00 


-1.9 5 


-1.66 




28.0 


-11.75 


-5.95 


-5.67 




27.0 


-15.44 


-9.87 


-9.61 
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as measured in the volumetric jar is taken as correct, then 
the per cent error in volume estimated by Model ffl can be in- 
dicated as a function of ot. The results are depicted in 
Figure (4.1). 

Before analyzing these results it would seem appro- 
priate again to describe the manner in which the measure- 
ments were taken. A rectangular grid boundary, 24 inches by 
60 inches was drawn around the base of the stockpile. Thus, 
sixty-one cross sections could be selected, at a one inch 
spacing, the first and last cross sections having zero cross 
sectional area. By measuring the distance between the base 
of the stockpile and the side boundary for each cross section, 
the base length of each individual cross section was deter- 
mined. The accuracy of these measurements was approximately 
l/lO of an inch because the scale was ruled in tenths and 
because individual particles were not much smaller than l/lO 
inch. 

From Figure (4.1a), it may be observed that the three 
quadrature rules give essentially the same error curve for 
volume determinations when the maximum number (61 ) of sec- 
tional areas is used in the computation. From the graph, 
the average apparent angle of repose for coated ammonium ni- 
trate granules was determined to be approximately 29.7 de- 
grees. The angle was taken to be the point at which the per 
cent error in the volume determination is zero. This appa- 
rent angle of repose would be used in Model #1 to estimate 
the volume of an actual undisturbed stockpile having 
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Figure 4.1 

Error in Volume Estimation by Model Number I 
for on Undisturbed Stockpile 
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essentially the shape of the one used in this particular 
test. The stockpile could consist of one or more peaks all 
essentially in line. 

Further study of Figure (4.1) reveals that the error 
curves for the three integration rules used in Model #1 di- 
verge as the number of cross sections is decreased. In order 
to study the effect of varying the number of cross section 
measurements on the volume calculation, additional volume cal- 
culations were made for 21, 16, 11, 6, 5, 4, and 3 cross sec- 
tions. These results together with the results for 61, 31, 

13, and 7 cross sections from Table (4.1) are listed in 
Table (4.2). 

Table (4.2) presents the per cent error in the calcu- 
lated volume utilizing the apparent angle of repose of 29.7 
degrees. The results are presented graphically in Figure (4.2). 
The per cent error in the volume behaves as a variable under- 
going damped oscillation. The error is large in magnitude 
for a small number of cross sections becoming damped to a 
small constant error as the number of cross sections ap- 
proaches one-half of the total available. Specifically, the 
figure indicates that the error curves converge at approxi- 
mately 31 cross sections and are quite close when 15 to 31 
cross sections are used. Thus, to compute the volume of a 
stockpile having essentially the same shape as the stockpile 
used in the present laboratory study one should use at least 
20 cross sections and possibly as many as 31 if the volume is 
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TABLE 4.2 



PER CENT ERROR IN CALCULATED VOLUME FOR AN 
APPARENT ANGLE OP REPOSE OF 29.7 DEGREES FOR 

MODEL PI 



Per Cent Error in Volume 


Cross Section/ 


Trapezoidal 


Simpson' s 


Weddle ’ s 


Interval Size* 


Rule 


Rule 


Rule 


61/1 


0.00 


0.00 


0.00 


31/2 


0.20 


-0.05 


-0.52 


21/3 


-0.82 


-0.21 


— 


16/4 


0.31 


— 


— 


13/5 


0.20 


2.10 


1.97 


11/6 


-2.52 


-6.04 


— 


7/10 


-5.30 


0.90 


1.20 


6/12 


6.60 


— 


— 


5/15 


-16.20 


-4.07 


— 


4/20 


-24.70 


— 


— 


3/30 


-52.20 


-36.80 


— 



* The notation 61/1, 31/2, etc., indicates 1 
interval containing 61 points, 2 intervals 
containing 31 points each, etc. Note that 
the terminal point on one interval is the 
initial point on the next interval, thus 
(points -1) X no. of intervals = 60 for all 
sets . 
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Figure 4.2 

Error in Volume Estimation by Model Number I 
for on Undisturbed Stockpile 
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to be estimated by Model #1. The appropriate apparent angle 
of repose would, of course, be 29.70 degrees. 

It is believed that the results of the above labora- 
tory model can be scaled to the dimensions of commercial 
stockpiles by considering the number of peaks in the commer- 
cial undisturbed stockpile. Thus, assuming the minimum num- 
ber of cross sections is 20 for the 4 peaked pile studied 
above, (see Figure 3.2) one might generalize the results to 
a rule that a minimum number of 5 cross sections per peak, 
i.e., 20 sections/4 peaks, should be used to estimate the 
volumes of undisturbed multi-peaked stockpiles of the type 
studied herein. 

It is to be noted that the Trapezoidal Rule gives a 
relatively low error level over the greatest range of cross 
sections. The Trapezoid Rule calculates the volume with an 
accuracy of - 1.0% for 13 or more cross sections while 
Simpson’s Rule calculates the volume to this same accuracy 
for 21 or more cross sections. Weddle’s Rule converges 
slower and therefore requires a greater number of cross sec- 
tions to compute the volume within - 1.07o. 

Model #2 

Table (4.3) tabulates the per cent error in the cal- 
culated volume obtained from computational Model #2 which 
involves two-way numerical integration. The following para- 
meters were varied in this study: 

1. number of cross sections 
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TABLE 4.3 

PER CENT ERROR IN CALCULATED VOLUME FOR UNDISTURBED 
STOCKPILE USING COMPUTATION MODEL #2 





Number of Depth Measurements 






25 PTS 


13 PTS 


7 PTS 


3 PTS 






% 


% 


% 


% 






Error 


Error 


Error 


Error 


Cross 


Integration 


Calc. 


Calc . 


Calc . 


Calc . 


Sections 


Rule* 


Volume 


Volume 


Volume 


Volume 




TRTR 


0.59 


0.60 


1.67 


2.57 




TRSN 


0.65 


0.67 


1.74 


2.65 




TRWD 


0.70 


0.71 


1.78 


2.72 


61 


SNTR 


0.59 


0.24 


2.35 


36.76 




SNSN 


0.65 


0.31 


2.41 


36.86 




SNWD 


0.69 


0.35 


2.46 


36.96 




• WDTR 


0.63 


0.16 


2.98 


99.94 




WDSN 


0.69 


0.23 


3.05 


100.08 




W)WD 


0.73 


0.27 


3.10 


100.22 




TRTR 


0.39 


0.39 


1.47 


2.34 




TRSN 


0.20 


0.24 


1.25 


2.14 




TRWD 


-0.23 


0.14 


0.75 


1.85 


31 


SNTR 


0.39 


0.02 


2.15 


36.45 




SNSN 


0.19 


0.10 


1.92 


36.18 




SNWD 


-0.26 


-0.44 


1.40 


35.80 




WDTR 


0.43 


-0.06 


2.78 


99.49 




WDSN 


0.22 


-0.15 


2.54 


99.02 




WDWD 


-0.24 


-0.45 


2.03 


98.36 




TRTR 


1.09 


1.11 


2.19 


2.65 




TRSN 


3.74 


3.82 


4.77 


6.13 




TRWD 


4.59 


4.70 


5.57 


7.62 


13 


SNTR 


1.08 


0.76 


2.53 


36.87 




SNSN 


3.72 


3.50 


5.20 


41.51 




S'NWD 


4.55 


4.41 


6.10 


43.49 




WDTR 


1.10 


0.66 


2.85 


100.13 




WDSN 


3.74 


3.42 


5.68 


106.59 




WDIVD 


4.58 


4.36 


6.74 


109.24 




TRTR 


-6.88 


-6.99 


-5.57 


-7.77 




TRSN 


0.64 


0.54 


2.02 


-0.33 




TRWD 


1.60 


1.52 


2.94 


0.88 




SNTR 


-6.85 


-7.46 


-5.48 


22.98 


7 


SNSN 


0.68 


0.05 


2.16 


32.90 




SNWD 


1.63 


1.04 


3.09 


34.51 




WDTR 


-6.82 


-7.62 


-5.63 


80.72 




WDSN 


0.70 


-0.10 


2.04 


95.17 




WDIVD 


1.63 


0.92 


3.01 


97.36 



*TR = Trapezoid Rule, SN = Simpson's Rule WD = Weddle's 
Rule, TRSN = Cross section integration using Trapezoid 
Rule and length integration using Simpson’s Rule 
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2. number of depth measurements per cross section 

3. quadrature rule for cross sections 

4. quadrature rule for length 

Because of the difficulty of conveniently represent- 
ing graphically all of the data contained in Table (4.3) only 
certain extreme examples are plotted herein. These corre- 
spond to the three curves which give consistently the small- 
est error in the volume calculation when using both the maxi- 
mum number of points (25 points) in the cross sectional area 
determination and the minimum number of cross section points 
(3 points). These results are plotted as Figure (4.3). 

It may be seen that the error curves have the same 

damped oscillation feature observed in Model with the 
error becoming relatively constant when approximately one 
half of the total possible cross sections on a one inch 
spacing are used, regardless of the number of depth measure- 
ments per cross section. As might be expected, the error is 
greater for the smaller number of depth measurements used per 
cross section. 

The results of Figure (4.3) show that a remarkably 
low error is incurred by using as few as three points per 
cross section when a large number of lengthwise cross sec- 
tions are taken. This error, less than 37o, is not appreci- 
ably greater than that obtained with a large number of cross 
sections and the assumption of an isosceles triangular cross 
section using Model # 1 , 
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Figure 4.3 



Error in Volume Estimation 
for an Undisturbed 
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This observation seems to indicate that, for Model ^2, 
one need not use a large number of data points in both the 
width and length dimensions on an undisturbed stockpile to 
obtain a satisfactory estimate of stockpile volume. Instead, 
one is probably justified in using a rather elementary ap- 
proach in obtaining cross sectional areas provided a sufficient 
number of these cross sections is used for the lengthwise in- 
tegration. 

An analysis of the calculated areas for two typical sec- 
tions was made to determine the best method of obtaining 
the area of any cross section. Table (4.4) tabulates the 
results of the area determinations for two typical cross sec- 
tions for the three quadrature rules. Using the maximum 
number of depth measurements, all quadrature rules computed 
the areas for the typical cross sections within .50%. 

Using the results obtained by applying Weddle's Rule 
to the maximum number of cross sectional depths measurements 
as the "correct” answer and then comparing the results ob- 
tained from both Simpson’s and the Trapezoidal quadrature 
Rules based upon only three depth raeasurmments , it may be 
observed that the Trapezoid Rule gives the lowest error. 

In a study similar to that shown in Table (4.4), but 
involving more cross sections, it was revealed that the 
errors in area estimation by the Trapezoidal Rule were oppo- 
site in sign for cross sections of high and of low area. 

These errors tended to cancel when used in the lengthwise 
integration formulas, thus further substantiating the use 



38 



TABLE 4.4 

COMPARISON OF COMPUTED CROSS SECTIONAL AREAS 
FOR UNDISTURBED STOCKPILE 



Parameters 


2 

Cross Sectional Area (in ) 


Number of Depth 
Measurements per 


Quadrature 


Number 13 


Number 7 


Cross Section 


Rule 


(large ) 


(small ) 


25 


Weddle ’ s 


74.66 


27.97 


25 


Simpson’ s 


74.62 


27.90 


25 


Trapezoidal 


74.62 


27.86 


3 


Trapezoidal 


72.00 


31.39 


3 


Simpson’ s 


96.00 


41.86 



Considering the 25- point Weddle’s Rule 
values to be the "correct” areas for 
each cross section the error at selected 
cross sections are as follows: 

Cross Section #13 

Per Cent Error in Area = 100.0 = -3.56% 

74.66 



. or Error in square inches = -2.66. 



Cross Section #7 

Per Cent Error in Area = 100.0 = +12.2% 



or Error in square inches 



+3.42. 
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of the simple Trapezoidal Rule for obtaining cross sectional 
areas. 

Table (4.5) summarizes for Model #2 the best alterna- 
tives in keeping with the basic purpose of this study - namely 
to determine the minimum number of measurement points and 
their location in order to best approximate the volume of a 
stockpile within the prescribed error range. The table pre- 
sents the alternatives for which the volume determination on 
the undisturbed stockpile, computed by Model #2, is less than 
- 3%. 

Prom the observations of the volume calculation methods 
for the undisturbed stockpile it may be stated that both 
Models #1 and #2 give excellent results when utilizing a 
maximum number of data measurements. Also, under appropriate 
circumstances, a much smaller number of measurements can be 
used in the models and still maintain the error levels within 
a tolerable range. Specifically, one is justified in using a 
rather simple geometrical representation in the cross-wise or 
cross sectional area determinations for elongated undisturbed 
stockpiles. It is possible to determine the cross sectional 
areas with a high degree of accuracy from as few as 3 to 5 
measurements, with approximately 5 cross sections being taken 
per peak over the length of the stockpile. 

If this procedure is followed using the Trapezoid or 
Simpson^ s quadrature Rule, then the stockpile volume should 
not be in error by more than - 3%. Further refinements are, 
of course, possible using additional measurements but 
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TABLE 4.5 

ALTERNATIVE MEASUREMENTS ON AN UNDISTORBED STOCKPILE 
FOR VOLUME ERROR RANGE OF - 3.0% 



Number of Depth 
Measurements per 
Cross Section 


Number 

of 

Cross Sections 


Total Number 
of 

Measurements 


Double 
Quadrature 
Ru 1 e * 


25 


61 


1525 


mm 

WDSN 

WDTR 


25 


51 


775 


WDSN 

WDTR 

WDWD 


25 


7 


175 


WDSN 

WDWD 


3(1) 


13 


39(11) 


TRTR 


3(1) 


7 


21(5) 


TRSN 

TRWD 



Note : 



Figures in parenthesis are actual 
number of measurements required 
since the end or outer boundary 
measurements are always zero. 



*The Code for the Double Quadrature Rule is as follows: 



WD = Weddle’s Rule 
SN = Simpson’s Rule 
TR = Trapezoidal Rule 

thus WDSN means Weddle’s Rule used in integrating 
the cross section and Simpson’s Rule used to inte- 
grate over the length. 
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accuracies greater than 1% are not to be expected. On the 
other hand, a very rapid approximation method using only 3 
points (in reality only 1 point) per cross section and 3 
cross sections per peak will give an error of less than 10%. 

B. Disturbed Stockpile Results 

Computational Models #1 and #2 were again used in 
obtaining the volume and apparent angle of repose of the 
disturbed stockpile. The method of analysis was essentially 
the same as that for Part A of this chapter. 

Model #1 

Table (4.6) presents the per cent error in the cal- 
culated volume obtained from computational Model #1 using 
the following as parameters; 

1. number of isosceles triangle cross sections 

2. quadrature rule 

3. apparent angle of repose. 

These results are shown graphically by Figure (4.4). It can 
again be observed that the error in the volume estimation by 
the three quadrature rules is essentially the same when the 
maximum number of cross sections is used in the computation. 
The graphs show that the apparent angle of repose for a dis- 
turbed stockpile is approximately 25.9 degrees. This value 
closely corresponds to the value of 26.5 degrees used by the 
Cooperative Farm Chemical Association in their calculations 
on disturbed stockpiles. Further observations of Figure (4.4) 
shows that the error curves are less widely separated for the 
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TABLE 4.6 

PER CENT ERROR IN CALCULATED VOLUME FOR 
DISTURBED STOCKPILE USING COMPUTATION 
MODEL ftl 





Per 


Cent Error in ' 


Calculated 


Volume 




Apparent 








Cross Sections 


Angle of 


Trapezoidal 


Simpson’ s 


Wedd le ’ s 




Repose 


Rule 


Rule 


Rule 




27.00 


4.80 


4.74 


4.41 




26.01 


— 


— 


0.00 


61 


26.0 


0.32 


0.26 


-0.06 


25.94 


— 


0.00 


— 




25.93 


0.00 


— 


— 




25.00 


-4.09 


-4.14 


-4.45 




27.00 


4.99 


5.31 


4.81 




26.00 


0.50 


0.81 


0.32 


31 


25.93 


— 


— 


0.00 


25.89 


0.00 


— 


— 




25.82 


— 


0.00 


— 




25.00 


-3.92 


-3.62 


-4.08 


■ 


27.00 


9.65 


10.07 


8.64 




26.00 


4.96 


5.36 


4.00 


13 


25.13 


— 


— 


0.00 


25.00 


0.35 


0.73 


-0.57 




24.92 


0.00 


— 


— 




24.84 




0.00 


““““ 


• 


27.00 


8.39 


16.29 


17.39 




26.00 


3.75 


11.32 


12.37 


f 


25.18 


0.00 


— 


— 




25.00 


-0.81 


6.43 


7.43 
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Figure 4.4 

Error in Volume Estimation by Model Number I 
for 0 Disturbed Stockpile 
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range of 7 to 61 cross sections than the corresponding curves 
applicable to the undisturbed stockpile. This implies a 
more nearly constant apparent angle of repose as the number 
of cross sections is decreased from the maximum of 61. This 
may be attributed to the increased homogeneity of the cross 
sections caused by the removal of material from the base. 

By reference to the disturbed stockpile. Figure (3.3) one 
can see that the peaks of the undisturbed stockpile have 
formed into ridges as a consequence of the removal of material 
from around the bottom of the stockpile. This action created 
a relatively planar stockpile surface in the region of the 
ridges. The forming of the ridges from each peak is there- 
fore, visual evidence of the homogeneous structure of dis- 
turbed stockpiles. 

Further volume calculations were made for the same 
additional cross sections as in Part A and the results are 
listed in Table (4.7) together with the applicable results 
of Table (4.6). Table (4.7) presents the per cent error in 
the calculated volume utilizing the apparent angle of repose 
of 25.9 degrees. 

Figure (4.5) depicts the error in the calculated 
volume as a function of the number of cross sections. For 
both the disturbed and the undisturbed stockpiles, approxi- 
mately 20 to 30 cross sections seem sufficient for an accu- 
rate volume determination with possibly slightly more cross 
sections being required for disturbed stockpiles since the 
curve for the disturbed stockpile does not seem to damp quite 
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TABLE 4-7 

PER CENT ERROR IN CALCULATED VOLu;U FOR 
APPARENT ANGLE OF REPOSE OF 25.9 DEGREES 



Cross Section/ 
Interval 


Trapezoidal 

Rule 


Simpson's | 

Rule 1 


Weddle'S ; 
Rule 


61/1 


.20 1 


.00 


-0.20 i 

j 


31/2 


.30 1 


.60 


0.10 ! 
1 


21/3 


-2.44 


-0.95 


-- i 


16/4 


-0.62 


— 


I 


13/5 


4.95 


5.35 


4.0 1 


11/6 


-6.80 


-5.20 


- j 


7/10 


3.75 


i 11.35 

I 


12.3 ’ 

1 


6/12 


-11.40 




' — i 


5/15 


-11.80 


-0.52 


i i 

! 


4/20 


-18.90 


— 


1 


3/30 


-69.00 


1 -28.4 


{ 

I i 
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Figure 4.5 

Error in Volume Estimotion by Model Number I 



for 0 Disturbed Stockpile 
(Angle of Repose = 25.9®) 
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as rapidly as for the undisturbed stockpile. It is not known 
however, how general is this difference in damping between 
the two types of stockpiles. For Model #1, an accuracy of 
- 3% is attainable with as few as 20 cross sectional measure- 
ments . 

In general, it may be concluded that the disturbed 
stockpiles require a slightly larger number of cross section 
intervals than do undisturbed stockpiles to compute the 
volume with the same accuracy. 

Model #2 

Table (4.8) presents the per cent error in the cal- 
culated volume obtained from computation Model #2 using the 
same parameters for Model #2 as used in Part A. A compari- 
son of these results with those listed in Table (4.3) for 
the undisturbed stockpile shows a general overall increase 
in the error of the computed volume for the disturbed stock- 
pile. This is especially noticeable when the number of 
depth measurements utilized in the computation, is decreased 
to a minimum. The error involving Weddle’s Rule using three 
depth measurements is an exception; the error value is mean- 
ingless with this small number of points since Weddle’s Rule 
requires 7 points. 

Figure (4.6) depicts the per cent error in the cal- 
culated volume obtained from the double quadrature rules 
which gave the minimum error calculation. Two sets of error 
curves are plotted; one utilizing the maximum number of 
depth measurements per cross section (23 points) and the 
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TABLE 4.8 

PER CENT ERROR IN CALCULATED VOLUME FOR DISTURBED 
STOCKPILE USING COMPUTATION MODEL #2 







Depth Measurements 






23 PTS 


13 PTS 


7 PTS 


3 PTS 






% 


% 


% 


% 






Error 


Error 


Error 


Error 


Cross 


Integration 


Calc . 


Calc. 


Calc . 


Calc . 


Sections 


Rule 


Volume 


Volume 


Volume 


Volume 




TRTR 


0.97 


1.25 


2.33 


9.53 




TRSN 


1.51 


1.80 


2.96 


9.97 




TRWD 


1.96 


2.26 


3.53 


10.37 




SNTR 


1.06 


0.89 


3.38 


46.05 


61 


SNSN 


1.61 


1.41 


4.04 


46.63 




SNWD 


2.06 


1.84 


4.63 


47.16 




WDTR 


1.09 


1.26 


6.86 


64.30 




WDSN 


1.64 


1.76 


7.53 


64.96 




WDWD 


2. 11 


2. 17 


8.13 


65.56 




TRTR 


0.31 


-0.39 


0.42 


9.65 




TRSN 


-1.10 


-0.79 


0.06 


9.20 




TRWD 


-1.77 


-1.47 


-0.61 


8.68 




SNTR 


-0.55 


-0.66 


1.39 


46.83 


31 


SNSN 


-0.15 


-1.07 


1.01 


46.61 




SNWD 


-1.88 


-1.75 


0.23 


46.22 




WDTR 


-0.72 


-0.22 


4.85 


64.75 




WDSN 


-1.14 


-0.66 


4.40 


64.22 




WDWD 


-1.88 


-0.34 


3.50 


63.36 




TRTR 


6.44 


1.87 


3.01 


9.65 




TRSN 


1.34 


3.01 


4.50 


11.10 




TRWD 


2.62 


3.09 


4.92 


12.52 




SNTR 


1.72 


1.49 


4.28 


46.90 


13 


SNSN 


2.53 


2.52 


6.02 


49.66 




SNWD 


2.46 


2.48 


6.67 


51.14 




WDTR 


1.54 


1.83 


7.92 


64.99 




WDSN 


2.51 


2.83 


10.03 


68.51 




WDWD 


2.41 


2.76 


10.94 


73.62 




TRTR 


-1.44 


-1.56 


-1.46 


5.00 


1 


TRSN 


3.65 


3.56 


3.73 


10.06 




TRWD 


3.22 


3.14 


3.31 


9.11 




SNTR 


-1.40 


-1.59 


-0.96 


40.00 


7 


SNSN 


3.69 


3.50 


4.26 


46.75 




SNWD 


3.25 


3.09 


3.68 


■ 45.48 




WDTR 


-1.35 


-1.17 


1.59 


57.50 




WDSN 


3.74 


3.93 


7.00 


65.09 




WDWD 


3.30 


3.54 


6.30 


63.37 
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Figure 4.6 

Error in Volume Estimation by Model Number 2 
for a disturbed Stockpile 
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other using the rainirauin number of depth measurements per 
cross section (3 points). The dampning characteristic for 
both sets of error curves is somewhat diminished compared to 
the corresponding curves (see Figure 4.3) for the undisturbed 
stockpile. The most noticeable trait of this plot is the 
relatively large difference between these two sets of error 
curves. The set of error curves using the maximum number of 
depth measurements per cross section shows a - 2% error for 
the range of 31 to 61 cross sections and an overall error 
level of - 6% for 7 to 61 cross sections. On the other hand, 
the set of error curves utilizing 3 depth measurements give a 
minimum error of approximately 10% over the complete range 
of cross sections. 

These results indicate that two-way numerical inte- 
gration over disturbed stockpile is not very accurate using 
a minimum number of depth measurements per cross section un- 
less one uses the Trapezoid Rule to determine the cross sec- 
tional areas and then integrate the resultant areas by either 
the Trapezoid or Simpson’s Rule to obtain the volume. 

It can therefore be concluded that for the disturbed 
stockpiles, which have relatively straight surfaces, the 
Trapezoidal quadrature Rule should be used for integration 
when using the minimum number of depth measurements per 
cross section. The volume error level which results using 
this method is reasonably stable but is of such a magnitude 
(10%) that another method for calculating the volume is to 
be preferred. 
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Table (4.9) presents for Model #2 for a disturbed 
stockpile, the best alternative measurements required to cal- 
culate the error in the volume to an accuracy of + 3.0%. A 
comparison of this table with the corresponding Table (4.5) 
in Part A shows that a greater number of depth measurements 
would be required on a disturbed stockpile to achieve the 
same 2 .% error level. 

This study has shown that numerical integration when 
using a minimum number of cross sectional depth measurements 
resulted in a error level that was 3 times greater than that 
determined by Model #1 using isosceles triangular cross sec- 
tions. It is therefore concluded that the use of the single 
isosceles triangle cross section is justified in approximat- 
ing the volume of disturbed stockpiles. 

In summary it may be stated that Model #1 seems pre- 
ferable for computing the volume of disturbed stockpiles. 

C. Angle of Repose Results 

The angle of repose of a granular stockpile is defined 
as the angle which a plane, tangent to the surface of the 
stockpile, makes with a horizontal plane. In this experiment, 
the undisturbed stockpiles consisted of sections that were 
conical and for this geometrical shape the angle of repose 
corresponds to the angle formed by the horizontal plane and 
by a line lying in the conical surface and passing through 
the apex of the cone. 

Since height measurements were taken at one inch in- 



tervals, in general, a cross section would not pass precisely 
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TABLE 4.9 

ALTERNATIVE MEASUREMENTS ON A DISTURBED STOCKPILE 
FOR A VOLUME ERROR RANGE OF - 3.0% 



Number of Depth 


Number 


Total Number 


Double 


Measurements per 


of 


of 


Quadrature 


Cross Section 


Cross Sections 


Measurements 


Rule* 








TRTR ' 


23 


61 


1403 


SNTR 

TRSN 








SNSN 


23 


31 


713 


TRTR 

TRSN 








TRSN 


23 


7 


161 

1 


SNTR 

SNIVD 

1 








TRSN 


7(5) 


31 


217(155) 


TRTR 

TRWD 


7(5) 


13 


91(65) 


TRTR 


7(5) 


7 


49(35) 


SNTR 

TRTR 



Note : 



Figures in parenthesis are actual 
number of measurements required. 
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through the apex of any of the four conically shaped peaks 
of the undisturbed stockpile. Therefore, in order to esti- 
mate the true angle of repose it was necessary to interpo- 
late between the cross sections bracketing the cross section 
which passed through the apex. 

A straight line was passed through the depth measure- 
ments, by the method of least squares, for each side of all 
cross sections. The slope of each line was determined and 
the angle less than 90 degrees corresponding to this slope 
was computed. These angles were then plotted versus the 
cross section position. The true angle of repose was then 
taken to be the interpolated angle corresponding to the 
known peaks cross section position. Figure (4.7) shows a 
typical plot of angles computed in the manner described above 
and corresponding to cross sections near as undisturbed stock- 
pile peak. The graph shows that the angle of repose is about 
31.6 degrees at cross section location 13.2. Table (3.1) 
lists the observed angle of repose by protractor measurement 
for cross section number 13 as 31.0 degrees. Figure (4.8) 
shows a plot of the height measurements for cross section 13 
together with the straight line obtained by the method of 
least squares. 

The angle of repose was also determined in the same 
manner, for the disturbed stockpile. Figure (4.9) is the 
plot of the angles of repose for the cross sections forming 
the first ridge section (observed at the left in Figure 3.3) 
Since the apex of the conical shape degenerated to a ridge. 
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Figure 4.7 

Angle of Repose for Undisturbed Stockpile 
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Figure 4.8 

Least Squares Fit of Depth Measurements 
for Cross Section Number 13 Undisturbed Stockpile 
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Figure 4.9 

Angle of Repose for Disturbed Stockpile 
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it can be seen that it is now possible to have adjacent angles 
of repose varying from 28.0 to 30.2 degrees for a disturbed 
stockpile, the average being approximately 29.5 degrees. 

One may think of the 31.6 degrees angle of repose for 
an undisturbed stockpile as the initial condition of a stock- 
pile. With the passage of time, factors such as weight of 
the granules, vibration, moisture, and size of granules cause 
the stockpile to compact and thereby reduce the angle of re- 
pose. For this experiment lasting approximately four weeks, 
the angle of repose decreased 2.1 degrees. 



CHAPTER V 



APPLICATION OF RESULTS 



General 

This chapter discusses the principles and methods of 
volume estimation for commercial scale stockpiles based on 
the methods developed in the previous chapters for laboratory 
scale stockpiles. The use of a digital computer program and 
a manual computation method utilizing the stockpile base 
measurements and the apparent angle of repose is explained. 

The computer program, Model #1, given in Appendix A, is 
adopted as the most appropriate mathematical model to be 
used in estimating the volume of stockpiles. This mathemat- 
ical model, which assumes that the typical stockpile cross 
section is an isosceles triangle, is the basis for the de- 
velopment of the manual computation forms shown in Appendix 
B. Appendix B contains one blank manual computation form 
and one completed sample calculation. 

Reserves Estimation Using a Digital Computer 

In using Model #1 computer program to compute the 
volumes of commercial scale stockpiles, it is first necessary 
to establish a boundary system in the stockpile warehouse. 

The boundary system can be either painted on the floor or 
associated with the stanchions of the warehouse. The side 
boundaries should be located as close to the base of the stock- 
pile as is physically possible. Also, the end boundaries 
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should be selected after proper recognition is taken of the 
number of cross sections to be used per peak. The total num- 
ber of cross sections taken should be an odd number. Further, 
although not absolutely necessary, it would be wise to select ' 
the interval between cross sections so that the cross sec- 
tions all fall on marked scale divisions of the base grid. 

These criteria allow a maximum flexibility in the use of the 
Trapezoid and Simpson’s integration formulas, and maximum 
simplicity in making the measurements. 

In determining the number of cross section base measure- 
ments to be taken, a decision needs to be made as to whether 
the stockpile is classified as "undisturbed” or "disturbed". 

If the stockpile is undisturbed, the general rule is to take 
at least three base measurements per each peak on the stock- 
pile. The resultant error in the estimated volume will be 
within approximately - 3%. If more than 3 measurements per 
peak are taken the resultant error in the estimated volume 
can be as low as - l7o. 

On the other hand, if the stockpile is disturbed, the 
general rule is to take at least 5 base measurements per each 
peak. The error in the estimated volume will be approximately 
in the range of - 3%. If more than 5 measurements per peak 
are taken the error in the estimated volume will approximate 
- 1 %. 

Before a final selection is made as to the number of 
measurements to be taken, consideration should be given to 
the rough estimate of the amount of material in the stockpile. 



60 



Good inventory management practice dictates that the larger 
the stockpile, the more accurate should be the volume deter- 
mination. Thus, the maximum number of base measurements 
should be taken consistent with the labor costs necessary to 
take the measurements. 

Once the number of base measurements to be taken has 
been determined, the next step is to begin measuring the 
distance between the base of the stockpile and the side ■ 
boundary line. The measurements should be to an accuracy of 
1/10 feet or the nearest inch. The cross section base measure- 
ments are required to be equally spaced over the length of 
the boundary grid and should of course be taken on both sides 
of the stockpile. The measurements should be recorded on a 
data sheet similar to that shown in Figure 1, Appendix B. 

Upon completion of the base measurements, the volume 
and tonnage estimates for the stockpile can be computed with 
a digital computer. The measurement data should be keypunched 
following the instructions given in Appendix A for the Model 
#1 computer program. The computer program also requires the 
length and width of the boundary grid together with the appa- 
rent angle of repose to be read in as data. The apparent 
angle of repose must be determined for different kinds of 
stockpile material in the manner outlined in Chapter III. The 
proper formats for all data cards are prescribed in the com- 
puter program. 

The computer program computes the estimated volume in 
the units selected for the input data. 
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The weight of material in the stockpile is then cal- 
culated using the average bulk density of the material. 

Reserves Estimation by Manual Calculation 

Appendix B contains both a blank manual computation 
form. Figure 1, and a completed manual computational form, 
Figure 2. This form was designed to use the same measure- 
ments as described for the computer program. The form is 
designed in three basic parts. At the top of the form is a 
block for constants. These are self explanatory and are 
used in the remaining two parts for computing the cross sec- 
tional areas, volume, and weight of stockpile material. The 
center section of the form calculates the areas of the cross 
sections utilizing the base measurement distances in Columns 
2 and 3 and assuming that a typical cross section is an isos- 
celes triangle. Columns 7, 8, and 9 of the form are totaled 
and these results form the basis for the volume and tonnage 
determination of the stockpile. 

In the computarion form, provision is made to compute 
the reserves by two different means; one using the Trapezoid 
Rule, the other using Simpson’s Rule for the lengthwise in- 
tegration of cross section areas. In general, these will 
give close but not identical answers as in the present case. 
It is suggested that either the two answers be averaged, or 
the most conservative estimate be selected and used for in- 
ventory management purposes. 



CHAPTER VI 



SUMMARY AND RECOMMENDATIONS 

In this thesis an attempt has been made to improve 
the volume estimation of commercial stockpiles through the 
judicious use of mensuration and numerical integration 
methods on laboratory scale stockpiles. 

During this study, the basic principle of mensura- 
tion was kept foremost in mind, namely, to determine the 
area of surfaces or the volume of solids from certain 
simple data on lines and angles. 

Using numerical analysis methods, it was found to be 
advantageous to unite the two computational methods. The 
experimental results show that the areas of all cross sec- 
tions of a stockpile can be readily and accurately approxi- 
mated by an isosceles triangle. Moreover, the area of an 
isosceles triangle can be determined by two simple measure- 
ments, its base length and angle of inclination. The appli- 
cation of either the Trapezoid or Simpson Rule of numerical 
integration to the approximate cross sectional areas re- 
sults in the most accurate volume estimates. The choice 
of the rule is somewhat arbitrary but depends essentially 
on the degree of regularity of the stock pile surface. It 
is recommended that the Trapezoid Rule be used for disturbed 
surfaces and the Simpson Rule for undisturbed surfaces, 
i.e., those formed by pouring material directly above the 
pile. Thus, advantage can be taken of both mensuration and 
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numerical integration concepts to satisfy the requirement for 
a minimum number of measurements to attain a relative high 
accuracy for the volume determination. 

The study has shown that, for a laboratory scale 
stockpile, this method of calculating the volume is realis- 
tic. It is believed that the results of the present labora- 
tory measurements can be adapted to the dimensions of com- 
mercial stockpiles by considering the number of cross section 
measurements taken per stockpile peak. For practical use, 5 
cross sections per peak are sufficient although 3 cross sec- 
tions per peak generally give errors of less than 10%. 

Finally, if this technique has awakened engineers and 
businessmen to the prospects of making a better determination 
of materials stored in a stockpile, for inventory management 
purposes, the main purpose of this thesis has been accomplished. 
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NOTES 



1. Homer A. Black and John E. Champion, Accounting in ^Busi - 
ness Decision , (Englewood Cliffs, N. J. , Prentice-Hall, 

Inc., 1962), p. 440. 

2. IBID . . p. 472. 

3. J. H. Perry, Chemical Engineers Handbook , 3rd Edition, 

(New York, McGraw-Hill Book Company, Inc., 1959) p. 1376. 

4. H. Stanley Redgrove, Experimental Mensuration , (London, 
England, Heinemann Publishing Co. LTD., 1912) Chapter XIII. 

5. J. B. Scarborough, Numerical Mathematical Analysis , 
(Baltimore, The Johns Hopkins Press, 1958), p. 131. 

6. IBID., pp. 51-127. 

7. IBID., p. 61. 

8. Leon Lapidus, Digital Computation for Chemical Engineers , 
(New York, McGraw-Hill Book Company, Inc., 1962), pp. 13-81. 

9. W. A. Granville, P. F. Smith and W. R. Longley, Elements 
of Calculus , (Boston, The Athenaeum Press, Ginn and Comp- 
any, 1946), pp. 113-167. 
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APPENDIX A 

COMPUTER PROGRAM TO DETERMINE THE APPRC..IMJITE VOLUME OF A 
STOCKPILE ASSUMING THE STOCKPILE CROSS SECTIONS ARE ISOSCELES 
TRIANGLES. 

Pur pose 

The purpose of this program is zo determine the approx- 
imate volume of a stockpile by numerical integration assuming 
isosceles triangle cross sections. A secondary purpose is to 
determine the angle of repose of the material forming the 
stockpile . 

LANGUAGE 

Fortran IV (IBM 7040 Computer) 

SYMBOLIC DICTIONARY 



VARIABLE 


S/A* 


I/O** 


DESCRIPTION 


XR 


A 


I 


Distance the base of the 
pile is from right bound- 
ary line. ^ 


XL 


A 


I 


Distance the base of the 
pile is from left bound- 
ary line . 


M 


S 


I/O 


Number of measurement 
points . 


L 


S 


I/O 


Distance betiveen the end 
boundary lines. 


W 


s 


I/O 


Distance between the side 
boundary lines. 



*S - Single variable; A - Array of variables 

■=^*1 - Input; 0 - Output 



I 

I 
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ALPHA 


S 


I/O 


Inc linr. cion angle of the 
stockpile - same as appar- 
ent angle of repose. 


ALPHAR 


S 


— 


Alpha expressed in radians. 


N 


S 


I 


Number of measurement 
stations . 


VACT 


s 


I/O 


Actual volume of the 
stockpile . 


DELTAL 


s 


0 


Length interval size. 


BASE 


s 


— 


Width of stockpile at the 
base . 


AREA 


A 


0 


Area of a cross section. 


SUM A 


s 




A holding variable used 
to accumulate the areas 
of cross sections for cal- 
culating the volume by the 
Trapezoidal Rule. 


VCALTR 


s 


0 


The calculated volume of 
the stockpile utilizing the 
Trapezoidal Rule. 


PERRTR 


s 


0 


Per cent error of the cal- 
culated volume using the 
Trapezoidal Rule. 


ODAREA 


s 




A holding variable used to 
accumulate the areas of the 
odd numbered stockpile 
cross sections for calculat- 
ing the stockpile volume us- 
ing Simpson’s Rule. 


EVAREA 


s 




A holding variable used to 
accumulate the areas of the 
even numbered cross sections 
for calculating the stock- 
pile volume using Simpson's 
Rule . 


VCALSN 


s 


0 


The calculated stockpile 
volume using Simpson's Rule. 


PERRSN 


s 


0 


Per cent error of the cal- 
culated volume using 
Simpson’s Rule. 
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SUMTWA 


S 




A holding variable used to 
accu.nulate the weighted 
areas of cross sections for 
calculating the stockpile 
volume using Weddle's Rule. 


VCALWD 


S 


0 


The calculated stockpile 
volume using Weddle's Rule. 


PRRRWD 


S 


O 


Per cent error of the cal- 
culated volume using Weddle's 
Rule . 



PROGRAM ROUTINE 

This program utilizes the data points (representing 
the distances to the side boundary lines) for each cross 
section. The area of each cross section is then computed 
using the angle of repose as a parameter. These areas are 
then integrated by the Trapezoid, Simpson's and Weddle's 
quadrature rules for a set interval size to obtain the cal- 
culated volume of the stockpile. The per cent error in the 
calculated volume is then computed. The angle of repose is 
determined when the per cent error is zero. 
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THIS PROGRAM DETERMINES THE APPROXIMATE VOLUME OF A STOCKPILE 
ASSUMING THAT THE CROSS SECTIONS ARE ISOCELES TRIANGLES. 

DIMENSION XL(IOO) » XR(lOO) » AREA(IOC) 

READ (5,5) M 

5 FORMAT ( 12 ) 

WRITE ( 6,6) 

6 FORMAT (15X, 20HM = NO. OF DATA PTS.//) 

WRITE (6,7) M 

7 FORMAT ( 19X, 12// ) 

READ (5,20) (XL(I) , XR ( I ) ,1 = 1, M) 

20 FORMAT (2F10.2) 

WRITE (6,11) 

11 FORMAT (15X, 17HINPUT DATA POINTS//) 

WRITE (6,13) 

13 FORMAT (15X, 3HNO., 15X, 2HXL, 15X, 2HXR//) 

WRITE (6,12) ( I» XL(I), XR(I) , I = 1, M) 

12 FORMAT ( 15X, 12, ICX, F10.2, 8X , F10.2//) 

9 READ (5,10) L, W, ALPHA, N, VACT 

10 FORMAT (2F10.2,F6.4, F6.0, F10.2) 

WRITE (6,14) 

14 FORMAT (15X, 20HINPUT DATA CONSTANTS//) 

WRITE (6,16) 

16 FORMAT(15X, IHL, lOX, IHW, 1 OX , 5HALPHA, lOX, IHN, lOX* 4hvAC' 
WRITE (6,17) L* W, ALPHA, N, VACT 

17 F0RMAT(13X, F5.2, 6X , F5.2, 8X , F5.2, 8X, F5.0- 9X, F10.2//) 

ASSUME CROSS SECTIONS ARE ISOCELES TRIANGLES FOR COMPUTING AREA 
real L,N 

ALPHAR = 3.1^159 ALPHA / 180.0 
C = ( SI N ( ALPHAR ) /COS (ALPHAR )) /4.0 

DELTAL = L/(N-1.0) 

WRITE (6,18) 

18 FORMAT (15X, 17HINTERVAL = DELTAL//) 

WRITE (6,19) DELTAL 

19 FORMAT (22X. F5.2//) 

J = N 

DO ICO I = 1, J 

BASE = W - XL ( I ) - XR ( I ) 

100 AREA(I) = C *BASF *BA5E 
WRITE (6,500) 

500 FORMAT (15X, 22HAREA OF CROSS SECTIONS//) 

WRITE (6,600) (I, AREA(I), I = 1, M) 

600 FORMAT(15X, 12* F10.4) 

STOCKPILE VOLUME COMPUTED WITH THE TRAPEZOID INTEGRATION RULE. 

K = J - 1 
SUMA = 0.0 
DO 80C I = 2» K 
800 SUMA = SUMA + AREA ( I ) 

VCALTR = (DELTAL/2.0) * (aREA(1) + 2.0 *SUMA + AREA(M)) 

PERRTR = 100.0 (VCALTR - VACT) / VACT 
WRITE ( 6,900) 

900 FORMATI15X, 5HALPHA, 15X- 4HVACT- 15X, 6HVCALTR, 15X, 6HPERRTR/ 
WRITE (6,400) ALPHA, VACT- VCALTR, PERRTR 
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400 FORMAT (15X, F5.2^ 15X- F7.2» 12X, F7.^, liX, Fb.2//) 

STOCKPILE VOLUME COMPUTED .V I THS IMPSONS INTEGRATION RULE. 

KK = J - 1 
2<K = KK - 1 
ODAREA = 0.0 
EVAREA = 0.0 
DO 910 I = 2. KK. ? 

910 EVAREA = EVAREA 4 - AREA (I) 

DO 920 I = 3» KKK» 2 
920 ODAREA = ODAREA + AREA ( I ) 

VCALSN = ( DELTAL/3 . ) * ( AREA { 1 ) + 4.*EVAREA + 2.*ODAREA + AREh(iM ) 
PERRSN = 100.0 (VCALSN - VACT) / VACT 
WRITE ( 6, 950) 

950 FORMAT (15X, 5HALPHA, 15X, AHVACTt 15X, bHVCALSN, 15X» 6HPERRSN// 
WRITE (6540^) ALPHA, VACT, VCALSN, PERRSN 

STOCKPILE VOLUME COMPUTED WITHWEDDLES INTEGRATION RULE. 

SUMTWA = 0.0 
DO 960 I ^ 2, M, 6 

960 SUMTWA = SUMTWA + l.*AR£A(I-l) + 5.*ARtA(I) + l.*AREA(I+l) + 

16. 4:-AREA ( 1+2 ) + 1.*AREA(I+3) + 5 . ’^-ARE A ( I +4 ! - 1 . -"-AR E A ( I + 5 ) 

VCALWD = . ?*DELTAL4;-SUMTWA 

PERRWD = 10U.0*(VCALWD - VACT) /VACT 

WRITE (6,970) 

970 FORMAT (15X, 5HALPHA, 15X, 4HVACT, 15X, 6HVCALWD, 15X, 6HPERRWL/. 
WRITE 16,400) ALPHA, VACT, VCALWD, PERRWD 
GO TO 9 
CALL exit 
END 
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FLOW diagram for COMPUTER PROGRAM TO DETERi^.INE THE APPROXIMATE 
VOLUME OF A STOCKPILE ASSUMING THE CROSS SECTIONS TO BE ISOS- 
CELES TRIANGLES. 





BASE = W - XL(I) - XR(I) 



AR£A(I) = C * BASE * BASE 
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PERRTR = 100.0*(VCALTR - VACT)/VACT 





0 



I 



75 





V 



a 



j 



76 




j/_ 



PERRSN = 100.0*(VCALSN - VACT)/VACT 



-V 



VVRITE 

ALPHA, VACT, VCALSN , PERRSN 



N/ 

VOLUME COMPUTED USING WEDDLE’S RULE 






SUMTWA = 0.0 
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COMPUTER PROGRAM TO DETERMINE THE APPROXIMATE VOLUME OF A 
STOCKPILE BY DOUBLE NUMERICAL INTEGRATION UTILIZING THE 
TRAPEZOID, SIMPSON’S AND WEDDLE’S QUADRATURE RULES. 

Purpose 

The purpose of this program is to determine the approxi 
mate volume of a stockpile by numerical integration of the 
height measurements for each cross section using the Trape- 
zoidal, Simpson’s and Weddle’s Rule. Then, using these re- 
sults, determine the volume of the stockpile through the use 
of the three quadrature formulas. Throughout the calculations 
the number of height measurements and the number of stockpile 
cross sections are designated as parameters. 

LANGUAGE 

Fortran IV (IBM 7040 Computer) 

SYMBOLIC DICTIONARY 



VARIABLE 


< 

CO 


I/O* ** 


DESCRIPTION 


CRSECT 


A 


I 


Height measurements in feet 
of the stockpile. 


AREATR 


A 


0 


Computed area of cross sections 
using the Trapezoid Rule. 


AREAS N 


A 


0 


Computed area of a cross section 
using Simpson’s Rule. 


AREAWD 


A 


0 


Computed area of a cross section 
using Weddle’s Rule. 

a 



*S - Single Variable; A - Array of Variables 

**I - Input; O - Output 



I 
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I 
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DIFF 1 


A 


0 


Difference in area computed by 
the Trapezoidal Rule and Simpson*; 
Rule . 


DIFP 2 


A 


0 


Difference in area computed by thi 
Trapezoidal Rule and Weddle's 
Rule . 


DIFF 3 


A 


0 


Difference in area computed by 
Simpson’s Rule and Weddle’s 
Rule . 


NPW 


S 


I 


Maximum number of cross section 
height measurements. 


NPL 


s 


I 


Maximum number of cross sectional 
areas for the stockpile. 


NPWX 


s 


I/O 


Actual number of cross section 
height measurements for a 
specific calculation. 


NPLX 


s 


I/O 


Actual number of cross sectional 
areas for a specific calculation. 


L 


s 


I/O 


Distance between the end 
boundary lines. 


W 


s 


I/O 


Distance between the side 
boundary lines. 


VACT 


s 


I/O 


Actual volume of the stockpile. 


DELTAW 


s 


0 


Cross section interval size. 


DELTAL 


s 


0 


Length interval size. 


SUMPTS 


s 




A holding variable used to 
accumulate the height measure- 
ments for calculating the areas 
of cross sections with both the 
Trapezoid Rule and Weddle’s Rule. 


ODDPTS 


s 




A holding variable used to accu- 
mulate the odd numbered . height 
measurements for calculating the 
areas of cross sections with 
Simpson’s Rule. 


EVNPTS 


s 




A holding variable used to accu- 
mulate the even numbered height 
measurements for calculating the 
areas of cross sections with 
Simpson’s Rule. 
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SUMA S 



ODAREA S 



EVAREA S 



VCTRTR S O 



VCTRSN S O 

VCTRWD S O 



A holding variable used to accu- 
mulate the areas of cross sections 
for calculating the volume of the 
stockpile utilizing the Trape- 
zoidal Rule and Weddle’s Rule. 

A holding variable used to accu- 
mulate the odd numbered areas of 
cross sections for calculating 
the volume of the stockpile 
utilizing Simpson’s Rule. 

A holding variable used to accu- 
mulate the even numbered areas 
of cross sections for calculating 
the volume of the stockpile 
utilizing Simpson’s Rule. 

The volume of the stockpile com- 
puted by numerical integration 
of the cross section height 
measurements by the Trapezoid 
Rule and then numerically inte- 
grating the resulting cross sec- 
tional areas by the Trapezoid 
Rule to obtain the volume of the 
stockpile . 

Volume computed using Trapezoid 
Rule and Simpson’s Rule. 

Volume computed using Trapezoid 
Rule and Weddle’s Rule. 



VCSNTR S O Volume computed using Simpson’s 

Rule and Trapezoid Rule. 

VCSNSN S 0 Volume computed using Simpson's 

Rule twice. 



VCSNWD S O Volume computed using Simpson’s 

Rule and Weddle’s Rule 



VCWDTR S O Volume computed using Weddle’s 

Rule and Trapezoid Rule. 

VCWDSN S O Volume computed using Weddle’s 

Rule and Simpson’s Rule. 



VCWDWD S 0 Volume computed using Weddle’s 

Rule twice. 
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PETRTR S 



PETRSN S 
PESNTR S 
PESNSN S 
PEWDTR S 
PEWDSN S 
PEWDWD S 



O The per cent error in the cal- 

culated volume when the calculated 
volume is obtained by numerically 
integrating the cross section 
height measurements with the 
Trapezoid Rule and subsequently 
numerically integrating the re- 
sulting cross sectional areas 
using the Trapezoid Rule. 

O The per cent error in the cal- 

culated volume utilizing Trape- 
zoid and Weddle’s Rules. 

0 The per cent error in the cal- 

culated volume utilizing Simp- 
son’s and Trapezoidal Rules. 

0 The per cent error in calculated 

volume utilizing Simpson’s Rule 
twice. 

O The per cent error in calculated 

volume utilizing Weddle’s and 
Trapezoidal Rules. 

O The per cent error in calculated 

volume utilizing Weddle’s and 
Simpson’s Rules. 

O The per cent error in calculated 

volume utilizing Weddle’s Rule 
twice . 



PROGRAM ROUTINE 

This program utilizes as data the height measurements, 
at one inch intervals, of the stockpile cross sections. The 
area of each cross section is then computed utilizing all 
three quadrature formulas while varying the number of height 
measurements from the maximum to the minimum. 

The resultant areas are then numerically integrated by 
three quadrature formulas to compute the volume of the stock- 
pile. The number of cross sectional areas is varied from the 
maximum to the minimum during this phase of the computations. 
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The resulting nine 
stockpile volume are then 
obtain the per cent error 



combinations for calculating 
compared to the known volume 
in the calculated volume . 
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THIS PROo^m.n CALCULaTL- THl V'Oi_UMii OF Ain. IRRLoOlaR STv^CKPIcE c ~^ 
OOU3LE IKTlGRATIOiN USING THE TRAPEZOIDAL. SIMPSONS. AND 
WEDDLES RULES. 



DIMENSION CRSECT ( 30.70 ) . AREATRI70). AREASNC’O). AREAw/D(70) 
DIMENSION DIFFK70). Dll-F2(70) . DIFF3I70) 

READ (5.5) NPW. NPL 

5 FORMAT (2131 
WRITE (6. A) 

A FORMAT ( m ) 

READ (5.26) NPWX , NPLX 
26 FORMAT (213) 

WRITE (6.6) 

6 FORMAT ( 1 5X ,22HNO. OF CRSECT DATA PT S . 5 X , 2 2HNO . OF LENGTH DATA PT 
WRITE (6.7) NPVvX . NPLX 

7 FORMAT (25X, 13. 2?X, 13//) 

READ(5.10) (( CRSECT ( I .J ) . 1=1, NPW). j-1. NPL) 

10 FORMAT ( 16F5.3/7F5.3) 



CONVERTING DATA VALUES INTO INCHES 
DO 1 1 J = 1 . NPL 
DO 11 I = 1. NPw 

11 CRSECT(I.J) = CRSECT! I ,J)-;U2.0 
READ (5.20) L. W. VACT 
20 FORMAT (213. F10.2 ) 



WRITE (6.21) 

21 FORMAT (15X, 20HINPUT DATA CONSTANTS//) 
WRITE (6.22) 

22 FORMAT (13X, IHL. 15X. IHW. 15X, AHVACT//) 
WRITE (6.23) L» W. VACT 

23 FORMAT (15X. 13. 12X. 13, 8X , F10.2//) 

INTFGFR W 

DELTAW = W/(NPWX - 1) 

DELTAL = L/(NPLX - 1) 

WRITE (6.2A) DELTAW, DELTAL 
2A FORMAT (15X, 9HDELTAW = . F5.2, 5X, 9HDELTAL 



F5.2// ) 



AREA OF CROSS SECTIONS COMPUTED BY TRAPEZOIDAL. SIMPSONS 
AND WEDDLES RULES 



1ST. TRAPEZOID RULE 



K = NPW - 1 
DO 30 J = 1 , NPL 
SUMPT5 = 0.0 
DO 31 I = 2 , K 

31 SUMPTS = SUMPTS + CRSECT! I, J) 

30 AREATR ( J ) = ( DELT AW / 2 . ) * ( CRSEC T ( 1 . J ) +2 . *S JMP TS+CRSEC T ( NPW . J ) ) 
2ND, SIMPSONS RULE 



KK = NPW - 1 
KKK = KK - 1 
DO 32 J = 1, NPL 
ODDPTS = 0.0 
EVNPTS = 0.0 

I 
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DO 33 I =2. KK, 2 

33 EVNPTS = EVNPTS + CRSECT(I-J) 

DO 34 M = 3, KK<» 2 

34 ODDPTS = ODDPTO + CRSECKM -J) 

32 AREASN(J) = ( DEL TAW/3 .)*( CRSECT ( 1 , J ) + 4.*EVNPTS + 2. -^ODDPTS + 
ICRSECT (NPaSJ) ) 

3RD, WEDDLES RULE 

DO 38 J = 1 , NPL 
SUMPTS = O.O 
DO 37 I =2, NPW, 6 

37 SUMPTS = SUMPTS + ] . -CRSEC T ( I - 1 , J ) + 5 . *CRSEC T ( I • J ) 

1+ 1 .*CRSECT ( I + l , J ) + 6.«-CRSECT( 1+2, J) + 1 . ^^'CRSECT ( I +3 , J ) + 
25.*CRSECT( 1+4, J 1 + 1 . *CRSECT ( I +5 , J ) 

38 AREAWD(J) = . 3«-SUMPTS-:^DELT AW 
DO ?'3 I = 1 , NPL 

DIFFl(I) = AREASN(I) - AREATR(I) 

DIFF2(I) = AREAWD(I) - AREATR(I) 

25 DlFF3(Ii = AREAWD(I) - AREASN(I) 

WRITE (6,35) 

35 FORMAT(15X, 6HAREATR* lOx, 6HAREASN, lOX, 6HAREAWD, lOX, 5HDIFF1 
15X, 5HDIFF2. 5X, 5HDIFF3//) 

WRITE (6,36) (I, AREATR(I), AREASN(I), AREAWD ( I ) ,D I FFl ( I ) , 

1CIFF2( I ) , DIFF3 ( I ) , I = 1, NPD 

36 FORMAT(10X, 13, FIO.^, 6X, F10.4, 6X, F1C.4, 5X , F7.4, 3X, FT. 4, 
13X, FT. 4) 

WRITE (6,39) 

39 FORMAT ( 3X ) 

COMPUTATION OF VOLUME USING TRAPEZOIDAL, SIMPSONS AND WEDDLES 
RULES. 

1ST, AREATR AND LENGTH USING TRAPEZOIDAL RULE 

N = NPL - 1 
SUM A = 0.0 

DO 40 I = ? , N 

40 SUMA = SUMA + AREATR(I) 

VCTRTR = (DELTAL/2.0) * (AREATR(l) + 2.*SUMA + AREATR(NPD) 
PETRTR = 100. 0*( VCTRTR - VACT)/VACT 

2ND, ARFATR AND LENGTH USING SIMPSONS RULE. 

N = NPL - 1 
NN = N - 1 
ODAREA = 0.0 
EVAREA = 0.0 
DO 50 I =2, N, 2 

50 EVAREA = EVARtA + AREATR(I) 

DO 51 I =3, NN, 2 

51 ODAREA = ODAREA + ARFATR(I) 

VCTRSN = ( DELTAL/3 .)*( AREATR ( 1 ) + 4. -^EVAREA + 2.*0DAREA + 

1 AREATR ( NPL ) ) 

PETRSN = 100.0*(VCTRSN - VACT)/VACT 
C 3RD, AREASN AND LENGTH USING TRAPEZOIDAL RULE 
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N = NPL - 1 
SUMA ^ 0.0 
DO 60 I = 2 » N 
60 SUMA = SUMA + AREASN(I) 

VCSNTR = (DELTAL/2. )-^(AREASN(l ) + 2.*SUMA + AREASN(NPD) 

PESiNTR = 100.0*(VCSNTR - VACT)/VACT 

4TH, AREASN AND LENGTH USING SIMPSONS RULE 

N = NPL - 1. 

NN = N - 1 
ODAREA = 0.0 
EVAREA = 0.0 
DO 70 I = 2» No 2 

70 EVAREA = EVAREA + AREASN(I) 

DO 71 I = 3o NNo 2 

71 ODAREA = ODAREA + AREASN(I) 

VCSNSN = ( DELTAL/3. )* (AREASN ( 1 ) + a.^^EVAREA + 2.*ODAREA + 
lAREASN ( NPL ) ) 

PESNSN = 100.0^>(VCSNSN - VACT)/VACT 

5THo AREATR AND LENGTH USING WEDDLES RULE 

SUMA = 0.0 

DO 75 I = 2 0 NPL 0 6 

75 SUMA = SUMA + 1 . *AR EAT R ( I - 1 ) + 5 . ^=-AREATR ( I ) + 1 . *AR EATR ( I + 1 ) + 

16. -x-aREATR ( I +2 ) (• 1 .">AREATR ( I +3 ) 5 . -AREATR ( I +4 ) + 1 . ->AR EAT R ( I +5 ) 

VCTRWD = . 3*SUMA*DELTAL 
PETRWD = 100.0* (VCTRWD - VACCj/VACT 

6TH, AREASN AND LENGTH USING WEDDLES RULE 

c,u«viA = 

DO 76 I = 2o NPLo 6 

76 SUMA = SUMA + 1 .* AREASN ( I - 1 ) + 5.*AREASN(I) + 1 . *ARFASN ( I + 1 ) + 

16. *AREASN ( 1+2 ) + i . *AREASN ( I +3 ) + 5 . *AREASN ( I +4 ) + 1 . *AR EASN ( I +5 ) 

VCSNWD = . 3*SUMA*DELTAL 

PESNWD = 100.0* (VCSNWD - VACT)/VACT 

7THo AREAWD AND LENGTH USING TRAPEZOIDAL RULE 

N = NPL - 1 
SUMA = 0.0 
DO 7 7 I = 2 0 N 

77 SUMA = SUMA + AREAWD(I) 

VCWDTR = ( DELTAL/2 .)*( AREAWD ( 1 ) + 2.*SUMA + AREAWD(NPL)) 

PEWDTR = 100.C*(VCWDTR -VACT)/VACT 

STHo AREAWD AND LENGTH USING SIMPSONS RULE 

N = NPL - 1 
NN = N- 1 
ODAREA = 0.0 
EVARE A = 0.0 
DO 78 I = 2o No 2 

78 EVAREA = EVAREA + AREAWD(I) 

DO 79 I = 3o NNo 2 
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79 ODAREA = ODAREA + AREAWO(I) 

VCWDSIM = (DELTAL/3. )«(AREAWD(1 ) + A.-*^EVAREA + 2.*ODAREA 
14- AREAWO(NPL) ) 

PEWDSi\ = 100.04;- (VCWDSN - VACT)/VACT 

9TH, AREAWD and LENGTH USING WEDDLES RULE 

SUMA = 0.0 

DO 30 I = 2» NPL» 6 

80 SUMA = SUMA 4- 1 . 44AREAWD ( I -1 ) 4- 5 . -»ARE AWD ( I ) t 1 . *ARE AWD ( I 4- 1 ) 4- 

16.-;4AREAWD( 1+2 ) +1 .-;4AREAWD( 1 + 3) + 5.*AREAWD( 1+4) + 1.44AREAWD( 1 + 5) 

VCWDWD = . 344SUMA44DELTAL 

PEWDWD = 100.0-)4(VCWDWD - VACT)/VACT 

WRITE (6, 41) VACT 

41 FORMAT (53X> lOHVACTUAL = ^ F10.2//) 

WRITE (6*42) 

42 F0RMAT(13X, 6 HVCTRTR> 6 X , 6 HVCTRSN» 6X , 6 HVCSNTR? 6 X » 6 HVCSNSN » 6 X 
1» 6 HVCTRWD> 6 X . 6 HVCSNWD 9 6 X • 6 HVCWDTR? 6 X» 6 HVCWDS\» 6 X> 6 HVCWDwD 
2» 6 X//) 

WRITE(6,52) VCTRTR» VCTRSN» VCSNTR* VCSNSN^ VCTRWD? VCSNWD» 
IVCWDTR? VCWDSN» VCWDWD 

52 FORMATdOX, Fl0.2» 2 X, F10.2? 2 X» F10.2» 2 X, F10.2> 2 X, F10.2» 2 X, 
1F10.2» 2 X, F10.2, 2 X, F10.2» 2 X> F10.2» 2 X//) 

WRITE (6»62) 

62 FORMAT (13X» 6 HPETRTR» 6 X. 6 HPETRSN» 6 X, 6 HPESNTR, 6 X, 

16HPESNSN» 6 X> 6 HPETRWD» 6 X » 6 HPESNWD* 6 X . 6 HPEWDTR. 6 X , 6 HPEWDSN* 
26X, 6 HPEWDWD. 6 X//) 

WRITE(6,72) PETRTR» PETRSN» PESNTR» PESNSId PETRWD* PESNWD, 
IPEWDTR? PEWDSNi PEWDWD 

72 F0RMAT(13X» F6.2. 6 X , F6.2» 6 X » F6.2, 6 X, F6.2? 6 X, F6.2» 6 X» 
1 F 6 . 2 ? 6 X, F6.2» 6 X , F6.2» 6 X , F6.2//) 

CALL EXIT 

END 



FLOW DIAGRAM FOR COMPUTER PROGRAM TO DETERMINE THE APPROXIMATE 
VOLUME OF A STOCKPILE BY DOUBLE NUMJERICAL INTEGRATION UTILIZING 
THE TRAPEZOIDAL, SIMPSON’S, AND WEDDLE’S QUADRATURE RULES. 
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XT 

M = 3 



<Sr 



ODDPTS = ODDPTS + CRSECT(M,J) 




AREASN(J) = (DELTAW/3.) * (CRSECT(1,J) + 

4. * EVNPTS + 2. * ODDPTS + CRSECTCNPW, J) 
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THIRD 

WEDDLE'S RULE 
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CALL EXIT 






END 
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Undisturbed Stockpile Depth Measurements 
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COMPUTER PROGRAM TO DETER.MINE THE APPROXIMATE VOLUME OF A 
BANKED STOCKPILE ASSUMING THE STOCKPILE CROSS SECTIONS ARE 
ISOSCELES TRIANGLES. 

Purpose 

The purpose of this program is to determine the approxi- 
mate volume o£ a stockpile} wrhoQa sides arc in contact with rc" 
taining walls. A secondary purpose is to determine the angle 
of repose of the material in the stockpile. 

Language 

FORTRAN IV (IBM 7040 COMPUTER) 



Symbolic Dictionary 



VARIABLE 


S/A* 


I/O** 


DESCRIPTION 


PW 




S 


I 


Distance between apex of the 
stockpile and the retaining wall. 


WF 




A 




Width of the small triangular 
section cut off by the retaining 
wall . 


XR 




A 


I 


Distance the base of the pile is 
from right boundary line. 


XL 




A 


I 


Distance the ba^e of the pile is 
from the left boundary line. 


M 




S 


I/O 


Number of measurement points. 


L 




S 


I/O 


Distance between the end boundary 
lines . 


W 




S 


I/O 


Distance between the side bound- 
ary lines . 


ALPHA 


s 


I/O 


Inclination angle of the stock- 
pile - same as apparent angle 
of repose. 


- 


Single 


variable 


; A - 


Array of variables 




Input ; 


0 - Output 
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ALPI-IAR 


S 


— 


Alpha expressed in radians. 


N 


S 


I 


Number of measurement stations. 


VACT 


S 


I/O 


. 0 

Actual’ volume of the stockpile. 


DELTAL 


s 


0 


Length interval size. 


BASE 


s 


— 


Width of stockpile at the base. 


AREA 


A 


0 


Area of a cross section. 


SUMA 


s 




A holding variable used to 
accumulate the areas of cross 
sections for calculating the 
volume by the Trapezoidal Rule. 


VCALTR 


s 


0 


The calculated volume of the 
stockpile utilizing the Trape- 
zoidal Rule. 


PERRTR 


s 


0 


Per cent error of the calculated 
volume using the Trapezoidal Rule 


ODAREA 


s 




A holding variable ’used to accu- 
mulate the areas of the odd num- 
bered stockpile cross sections 
for calculating the stockpile 
volume using Simpson's Rule. 


EVAREA 


s 




A holding variable used to accu- 
mulate the areas of the even 
numbered cross sections for cal- 
culating the stockpile volume 
using Simpson's Rule. 


VCALSN 


s 


0 


The calculated stockpile volume 
using Simpson’s Rule. 


PERRSN 


s 


0 


Per cent error of the calculated 
volume using Simpson’s Rule. 


SUMTIVA 


s 




A holding variable used to accu- 
mulate the weighted areas of 
cross sections for calculating 
the stockpile volume using Weddle 
Rule . 


VCALV® 


s 


0 


The calculated stockpile volume 
using Vu'eddle’s Rule. 


PERRV7D 


s 


0 


Per cent error of the calculated 
volume using 'leddle’s Rule. 



Ill 



Program Routine 

The program utilizes as data points the distances 
which the stockpile cross sections are from the side boundary 
lines. The area of each cross section is then computed assum- 
ing that it is an isosceles triangle. The program checks to 
see if the stockpile is in contact with the retaining wall. 

If it is, then the program calculates the net area of the 
cross section by subtracting the area of the cut-off tri- 
angular section from the gross area of the standard isosceles 
triangles . 

These areas are then integrated using the Trapezoidal, 
Simpson's and Weddle's quadrature formulas for a set inter- 
val size. The error in the calculated volume is then com- 
puted. 

The angle of repose is determined when the calculated 
volume error is a zero value based upon varying the parameter, 
ALPHA. 



0 
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THIS PROGRAM DETERMINES THE APPROXIMATE VOLUME OF A STOCKPILE 
WHICH IS BANKED AGAINST A RETAINING WALL AND FURTHER 
ASSUMING THAT THE CROSS SECTIONS ARE ISOSCELES TRIANGLES. 

DIMENSION XL(IOO), XR(IOO). AREA(IOO), BASE02(100)f WF(IOO) 
READ(5^26) PW 
26 FORMAT (F5.1 ) 

READ (5.5) M 

5 FORMAT ( 12 ) 

WRITE ( 6>6) 

6 FORMAT (15X, 20HM = NO. OF DATA PTS.//) 

WRITE (6»7) M 

7 FORMAT (19X, 12//) 

READ (5»20) (XL(I) > XR ( I ) . I = 1. M) 

20 FORMAT (2F10.2) 

WRITE (6.11) 

11 FORMAT (15X, 17HINPUT DATA POINTS//) 

WRITE (6.13) 

13 FORMAT (15X, 3HN0.. 15X, 2HXL. 15X. 2HXR//) 

WRITE (6.12) ( I. XL(I). XR(I) . 1 = 1, M) 

12 FORMAT ( 15X, 12. lUX. F10.2. 8X . F10.2//) 

9 READ (5.10) L. W. ALPHA, N. VACT 

10 FORMAT (2F10.2.F6.4. F6.0, F10.2) 

WRITE (6,14) 

14 FORMAT (15X, 20HINPUT DATA CONSTANTS//) 

WRITE (6.16) 

16 F0RMAT(15X, IHL. lOX. IHW. lOX. 5HALPHA, lOX. IHN, lOX, 4HVACT//) 
WRITE (6.17) L. W, ALPHA, N. VACT 

17 F0RMAT(13X, F5.2, 6X, F5.2. 8X , F5.2. 8X, F5.0, 9X, F10.2//) 

ASSUME CROSS SECTIONS ARE ISOSCELES TRIANGLES FOR COMPUTING AREAS. 

real l . n 

ALPHAR = 3.14159 ALPHA / 180.0 
A = SIN(ALPHAR)/C0S(ALPHAR) 

DELTAL = L/(N-1.0) 

WRITE (6,18) 

18 FORMAT (15X, 17HINTERVAL = DELTAL//) 

WRITE (6,19) DELTAL 

19 FORMAT (22X. F5.2//) 

J = N 

DO 44 I = 1 , J 
IF (XR( I ) .LE.0.0) GO TO 33 
22 BASE02(I) = ( W-XR ( I ) -X L ( I ) ) / 2 . 0 
AREA(I) = (BASE02( I )-^*2 )*A 
GO TO 44 

33 BASE02(I) = ( W-2 . O^f XL ( I ) ) /2.0 
WF ( I ) = BASE02( I ) -PW 

100 AREA(I) = ( BASE02 ( I ) **2 ) *A - ( WF ( I ) ^s-*2 ) *A*0 . 5 
J = I 

44 CONTINUE 

AREA! 1 ) = 0.0 
AREA(J) = 0.0 
WRITE (6,500) 

500 FORMAT (15X, 22HAREA OF CROSS SECTIONS//) 



n n n n 
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WRITE (6»600) (I, AREA(I), I = 1, M) 

60C F0RMAT(15X, 12, FIO.A) 

C STOCKPILE VOLUME COMPUTED WITH THE TRAPEZOID INTEGRATION RULE. 

K = J - 1 
SUMA = 0.0 
DO 800 I = 2, K 
800 SUMA = SUMA + AREA(I) 

VCALTR = (DELTAL/2.0) * (AREA(l) + 2.0 *SUMA + AREA(M)) 

PERRTR = 100.0 * (VCALTR - VACT) / VACT 
WRITE ( 6,900) 

900 F0RMAT(15X, 5HALPHA, 15X, 4HVACT, 15X, 6HVCALTR, 15X, 6HPERRTR//) 
WRITE (6,400) ALPHA, VACT, VCALTR, PERRTR 
400 FORMAT (15X, F5.2, 15X, F7.2, 12X, F7.2, 13X, F6.2//) 

STOCKPILE VOLUME COMPUTED WITHSIMPSONS INTEGRATION RULE. 

KK = J - 1 
KK'< = KK - 1 

ODAREA = 0.0 
EVAREA = 0.0 
DO 910 I = 2, KK, 2 
910 EVAREA = EVAREA + AREA (I) 

DO 920 I = 3, KKK, 2 
920 ODAREA = ODAREA + AREA ( I ) 

VCALSN = (DELTAL/3. )*(AREA( 1 ) + 4.*EVAREA + 2.*ODAREA + AREA(M)) 
PERRSN = lOo.O * (VCALSN - VACT) / VACT 
WRITE ( 6, 950) 

950 FORMAT(15X, 5HALPHA, 15X, 4HVACT, 15X, 6HVCALSN, 1 5X , 6HPERRSN//) 
WRITE (6,400) ALPHA, VACT, VCALSN, PERRSN 

STOCKPILE VOLUME COMPUTED WITHWEDDLES INTEGRATION RULE. 

SUMTWA = 0.0 
DO 960 I = 2, M, 6 

960 SUMTWA = SUMTWA + l.*AREA(I-l) + 5.*AREA(I) + l.*AREA(I+l) + 

16. *AREA ( 1+2 ) + 1 .^^AREA ( 1+3 ) + 5.*AREA(I+4) + 1 . *AREA ( I +5 ) 

VCALWD = .3*DELTAL-«SUMTWA 
PERRWD = 100.0*(VCALWD - VACT) /VACT 
WRITE (6,970) 

970 FORMAT(15X, 5HALPHA, 15X, 4HVACT, 15X, 6HVCALWD, 15X, 6HPERRWD//) 
WRITE (6,400) ALPHA, VACT, VCALWD, PERRWD 
GO TO 9 
CALL EXIT 
END 
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FLOW DIAGRAM FOR COMPUTER PROGRAM TO DETERMINE THE APPROXIMATE 
VOLUMJ: OF A BANKED STOCKPILE ASSUMING THE CROSS SECTIONS TO BE 
ISOSCELES TRIANGLES. 





READ 

PW 



X'-' 



READ 

M 



READ 

XL(I), XR(I), 
I = 1, M 



READ 

L, W, ALPPIA, N, VACT 



REAL 


L, N 







ALPHAR = 3.14159 * ALPHA/180.0 



A = SIN (ALPHAR )/COS (ALPHA R) 
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FIGURE I 

CALCULATION OP STOCKPILE VOLUME FROM BASE MEASUREr^NTS 



Boundary Width 


Anr;le of Rcoose 


Tan oC » 


C ^ Tan oC. /4,0 


Int crvp] Slee^I 


Boun-i'iry lenrjth 


RV = ft. 


- 




C « 


«BL/()» of CrosfJ 
'^ect , -1 ) = 


BI » ft. 



Col.(l) 


Col. (2) 


Col. (3) 


Col. W 


Col. (5) 


Col. (6) 


Col. (7) 




Co1.(r) 


Col. (9) 


Cross 

Section 


Distance to 
Pile Base 
from Ht, Sid 
Ft. 


Distance to 
Pile Base 
i from Lt.Side 
Ft. 


.Co1(2)4-Co1(: 

Ft. 


) Width of 
Pile Base 
BW -Col W 
Ft. 


Col (5) X 
Col (5) 

Ft? 


Cross .Sect. 
Area 

Col (6) 3CC 
FtT 




:ODD CiCvr/:?) 
X 2 

Col (7) ODD 

X 2 


fr.V vN Ca:t^ 

X ^ 

Col (?) /ViJJ 
X U 






































































































































































































0 






































































































































































































































































TOTAL 





















ESTIMATE USING TRAPEZOID RULE 


ESTHETE USING SIMPSON’S RULE 


Volume • Interval x Total Col (?)« 


Volume - Interval x Total Col (B)-*- Total Col (Q) ® 


Tonnaffe » Volume x Bulk Density X ( 1 1 « 

(2000 ) 


Tonnpee * Volume x Bulk Density x ( 1 ) - 

(2000 ) 
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FIGURE 2 

CALCUIATION or STOCKI’II.K VOl.UMK FUOM HASi: MKASUUUMKNTS 
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TOTAL 












854.38 
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ESTIMATE USING TRAPEZOID RULE 


ESTIMATE USING SIMPSON’S Rllld-: 


Volume* Interval x T 0 UI Co1(7)-(1*0) (854.38) » 854.38 C0.F1 


Volume* Interval xCTotal Col ( 8 ) ‘♦‘Total Col (9)) 
3 

3 U) (852.58 *’' 1712.36) - 854.98 CU.FT. 

ni 


Tftnnj»c*es Volume X Bulk Density! ( 1 

(2000 ) 

• (95'^.38)('<7a9)(_l_) 3 J0.16 T0^'.‘; 
(nnnn 1 


TonnsAe = Volume x Bulk Den-ity x (_ 1 1 

(2000 ) 

• (854.98)(47.19)( 1 ) „ 20.16 TONS 

(2000 ) 
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